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Abstract — In contrast to the classical cyclic prefix (CP)-OFDM, 
the time domain synchronous (TDS)-OFDM employs a known 
pseudo noise (PN) sequence as guard interval (GI). Conventional 
channel estimation methods for TDS-OFDM are based on the 
exploitation of the PN sequence and consequently suffer from 
intersymbol interference (ISI). This paper proposes a novel data- 
aided channel estimation method which combines the channel 
estimates obtained from the PN sequence and, most importantly, 
additional channel estimates extracted from OFDM data sym- 
bols. Data-aided channel estimation is carried out using the 
rebuilt OFDM data symbols as virtual training sequences. In 
contrast to the classical turbo channel estimation, interleaving 
and decoding functions are not included in the feedback loop 
when rebuilding OFDM data symbols thereby reducing the 
complexity. Several improved techniques are proposed to refine 
the data-aided channel estimates, namely one-dimensional (1- 
D)/two-dimensional (2-D) moving average and Wiener filtering. 
Finally, the MMSE criteria is used to obtain the best combination 
results and an iterative process is proposed to progressively 
refine the estimation. Both MSE and BER simulations using 
specifications of the DTMB system are carried out to prove 
the effectiveness of the proposed algorithm even in very harsh 
channel conditions such as in the single frequency network (SFN) 
case. 

Index Terms — TV broadcasting, channel estimation, OFDM, 
iterative method. 



I. Introduction 

CONVENTIONALLY, orthogonal frequency division 
multiplexing (OFDM) adopts cyclic prefix (CP) as guard 
interval (GI) to mitigate intersymbol interference (ISI) and 
enable simple equalization at receiver side. However, CP is 
a redundant copy of the data, which reduces the useful bit- 
rate and degrades the spectrum efficiency. Recently, it has 
been proposed to replace CP by a known pseudo-noise (PN) 
sequence. The main advantage of this approach is that the 
known sequence, which initially serves as GI, can be reused 
as training sequence to carry out channel estimation and syn- 
chronization. Therefore no pilot is inserted among the useful 
data symbols in the frequency domain as needed in traditional 
CP-OFDM. This kind of OFDM waveform is referred to as 
time domain synchronous-OFDM (TDS-OFDM£| Q and has 
been adopted by the Chinese digital television/terrestrial mul- 
timedia broadcasting (DTMB) system [4|, [5]. After perfectly 
removing the PN sequence, the received TDS-OFDM signal 
is converted to the so-called zero padding (ZP)-OFDM |6|. 
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'it is also known as pseudo random postfix-OFDM (PRP-OFDM) [2] and 
known symbols padding-OFDM (KSP-OFDM) | 3 | in other literatures. 



Then, all the well established equalization techniques 10, Q 
for ZP-OFDM can be applied for TDS-OFDM signals. 

The usage of the PN as GI is somehow "double-edged 
sword". On one hand, it improves the specttum efficiency by 
economizing the cost of transmitting redundant data in the 
CP and the pilots in the data symbols; on the other hand, 
the channel memory effect introduces mutual interference 
between the PN sequence and the OFDM data symbols. That 
is to say, the interference from data symbols compromises the 
performance of the PN-sequence-based channel estimation and 
the inaccurate channel estimate consequently causes imperfect 
PN sequence removal and further increases the difficulty to 
recover the data symbols. The channel estimation is thus the 
key point in the TDS-OFDM based transmissions. 

Many papers have addressed this challenging issue in the 
past few years and can be classified into three categories: 

1. Making channel estimation based only on the PN se- 
quence. In Q3, the time domain cross-correlation between the 
received PN sequence and the locally generated one is used to 
make channel estimation. It exploits the correlation property of 
the PN sequence in order to obtain processing gain to combat 
the noise and ISI. The complexity significantly grows with the 
increase of the length of the channel delay spread. In 0, a 
simple channel estimator is proposed based on the expectation 
(mean value) of the received PN sequence which cancels the 
random ISI component. However, this algorithm only adapts 
to the static channel situations and requires large storage. 

2. Removing ISI from the PN sequence when making channel 
estimation. A first approach [ 1 ] proposed to iteratively remove 
the ISI using the decision feedback method. However, much 
complexity is spent by a series of fast Fourier transforms 
(FFT). In |9|, it is proposed to carry out the channel estimation 
using the PN sequence and its tail spread in the beginning part 
of each OFDM data symbol. By iteratively removing the data 
symbol, the channel estimation is progressively improved. The 
estimation results however suffer an error floor even at low 
signal to noise ratio (SNR) due to the existence of ISI in the 
PN sequence. JT0| proposed to alternately remove the mutual 
interference between the PN sequence and the OFDM data 
symbols. A so-called partial decision is also used to reduce 
the error detection of data symbols. However, the performance 
significantly degrades when the channel delay spread is long. 
ifTTI also adopted the partial decision technique to remove the 
interference on the PN sequence from previous data symbols. 
In order to avoid the interference from the following data 
symbols, the "tail" of the received PN sequence is estimated 
and then used to reconstruct the interference-free received PN 
sequence for channel estimation. 
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Fig. 1. Block diagram of TDS-OFDM system. The shaded blocks are the processing used for the proposed method. The dashed blocks are the additional 
processing needed for the turbo channel estimation. 



3. Using the pilots in a complementary way to the PN 
sequence to enhance the channel estimation. It is proposed 
to replace some data subcarriers by pilots to improve the 
channel estimation performance in [3|. However, the spectrum 
efficiency is compromised. Another possible solution is the 
data-aided channel estimation method which uses the rebuilt 
data symbols as "virtual pilots" to make channel estimation 
and thus does not need any extra pilots. In lfl2ll . it is proposed 
to make hard decisions to the equalized symbols that fall in the 
reliable decision region. The hard decided data symbols are 
used as virtual pilots for channel estimation. However, this 
method requires an additional PN frame header for several 
OFDM symbol to obtain a good initial channel estimate. 
Therefore, it is not straightforwardly applicable to the existing 
DTMB system. In [13|, a turbo channel estimation algorithm 
is proposed in the context of CP-OFDM. An obvious disad- 
vantage of this method is the extremely high complexity and 
long time delay, especially in the systems with sophisticated 
channel decoder and deep interleaves For example, the DTMB 
system adopts LDPC codes and convolutional interleaver with 
a time delay of 170 (or even 510) OFDM frames which 
prohibits the using of such a method. 

The aim of this paper is to propose a low-complex but 
effective algorithm to increase the robustness of the channel 
estimation in TDS-OFDM compared to the already existing 
techniques described above. The highlights of algorithm pro- 
posed in this paper are (a) excluding the channel decoder, 
interleaver and de-interleaver from the iterative channel es- 
timation process to reduce complexity, (b) rebuilding soft 
symbols for channel estimation to prevent error propagation, 
(c) exploiting the correlation property of the channel to obtain 
improved estimation results, and (d) minimum mean square 
error (MMSE) combination of channel estimates obtained 
from PN and data. 

More concretely, we propose to rebuild the data symbols 
using the likelihood information from the demapper. Then the 
instantaneous data-aided channel estimates are acquired using 
these rebuilt data symbols. As the interleaving and channel 
decoding processes are not included in the feedback loop, the 
computational complexity is significantly reduced compared to 
the turbo channel estimation. The soft symbols are fed back in 
the iterations serving as virtual pilots in the channel estimation. 
Since no decision is made on the rebuilt symbols, the proposed 



method does not cause error propagation due to imprudent 
hard decisions. The uncertainty (noise) is kept in the soft 
data symbols and in the resulting data-aided channel estimates 
as well. Moreover, we propose several improved techniques, 
namely one-dimensional (l-D)Awo-dimensional (2-D) moving 
average and Wiener filtering, which exploit the time-frequency 
correlation property of the channel to suppress the noise ex- 
isting in the data-aided channel estimates. The cooperation of 
these refining techniques and soft symbol rebuilding achieves 
satisfactory channel estimation accuracy with low computa- 
tional complexity. The final channel estimates combine the 
results obtained from PN and from data according to MMSE 
criteria. The proposed method can adapt to difficult situations, 
including higher order constellations, extremely long channel 
spread and channel time variations. 

The rest of the paper is organized as follows. In section 
II, the mobile channel and TDS-OFDM signal models are 
described. In section III, the PN based channel estimation is 
presented. The new proposed data-aided channel estimation 
method is presented in section IV. Computational complexity 
is evaluated in section V. Simulation results are shown in 
section VI. Conclusions are drawn in section VII. 

In this paper, (•)* denotes the conjugate of complex number, 
E{-} is the expected value, (-) T and (-) H are the matrix 
transpose and Hermitian transpose, respectively. 

II. System Model 

A. Discrete Channel Model for OFDM 

The wireless channel is modeled as an L th order time- 
varying finite impulse response (FIR) filter. The parameter L 
is determined by the maximum excess delay of the channel. 
The channel is assumed to be quasi-static, namely channel co- 
efficients remain constant within one OFDM symbol duration 
but change from one OFDM symbol to another. The channel 
impulse response (CIR) for the i th OFDM symbol is: 

L-l 

h[i,m] = ^2h^S[m-l], (1) 

1=0 

where 5[-] is the Rronecker delta function, of 1 is the I th filter 
tap and is modeled as zero mean complex Gaussian random 
variables. Furthermore, the channel filter taps are assumed 
to follow the wide sense stationary uncorrelated scattering 
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Fig. 2. PN removal of the TDS-OFDM signal. 



(WSSUS) assumption, namely different taps are statistically 
independent, while a specific tap is correlated in time. The 
channel frequency response (CFR) for the fc th subcarrier of 
the i lh OFDM symbol can be obtained via an TV-point FFT 
over the CIR: 



H[i,k] 



N-l 

E 

1=0 



Ilk 



(2) 



The time-frequency correlation function of the CFR is lfT31 : 

</>H [P, q] = E {H[i + p, k + q]H* [i,k)} = r t [p]r f [q] , (3) 

where r t [p] and r/[g] are the time and frequency domain 
correlation functions of the CFR, respectively. Concretely, 
taking the Jakes' mobile channel model, r t [p] is [15|: 



n\p] = Jo(2Kpf d T b ), 



(4) 



where Jo(-) is the zero-order Bessel function of the first kind, 
and fd is the maximum Doppler frequency related to velocity 
v and carrier frequency f c by fd = vf c /c, where c is the speed 
of light, Tf, is the time duration of one OFDM block and T), = 
T g + T given T g and T the durations of GI and OFDM data 
parts, respectively. The frequency domain correlation rf[q] is: 



L-l 



r f[q\ 



(5) 



1=0 



where of is the power of the I th path. Without loss of 
generality, the power of the channel is normalized so that 

EL-1 2 i 
i=o «T = !• 

B. TDS-OFDM Signal Model 

Fig. Q] presents the baseband model of the TDS-OFDM 
based system. The i th OFDM data symbol is formed by TV- 
point inverse fast Fourier transform (IFFT): 
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A j/-length PN sequence {cP' [ri\}n=-v * s trien inserted 
before {a;^^^]},^) 1 as GI. The i th transmitted time domain 
signal is thus: 
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Fig. 3. Illustration of the OLA process. 



After passing the multipath fading channel, the received 
signal is a linear convolution of the transmitted signal and 
the channel as given in (0 which is shown in the top of next 
page, 

where [n]jv is the residue of n modulo N and w is the 
additive white Gaussian noise (AWGN) with a variance of 
crj. From (O and Fig. [2] it can be found that there are some 
mutual interference between the OFDM data symbols and the 
PN sequences due to the channel memory. In order to recover 
the data symbols, it is necessary to remove the PN sequence 
from the received signal first. Since the PN sequence [n] 
is perfectly known by the receiver, the channel-distorted PN 
sequence can be estimated by making a linear convolution of 
PN sequence c^[n] and the estimated channel response h^: 
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-v < n < L — 1. 



(9) 



Using cW[n],the PN sequence is removed from the received 
signal as depicted in Fig. If the channel is perfectly 
estimated, the PN sequence and its tail can be completely 
removed from the received signal. Otherwise, there will be 
some residual contribution of the PN sequence in the received 
signal. The PN sequence removal is expressed in (TTOb which 
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is given in the next page, where Ah l 
estimation error of the I th channel tap. 

After removing the PN sequence, the TDS-OFDM signal is 
turned to an equivalent ZP-OFDM signal. Thus, the so-called 
overlap and add (OLA) [3| process is performed by adding 
the following GI to the beginning part of the OFDM symbol 
as shown in Fig. [3] More precisely, the OLA process can be 
written as: 

y (' i; )(n) = f«(n)+f (l+1) (n-;y), < n< L - 1. (11) 

Suppose that the length of GI exceeds that of the channel, 
there is no ISI between two adjacent OFDM symbols. Taking 
( [Tol l into account, the received signal after removing the GI is 
written as: 



L-l 



°N=E x (l) [n - l] N + w'[n]+(,[n],0 < n < N, (12) 



1=0 



where w' (n) is the noise after OLA which is slightly colored 
and boosted in the OLA process. The equivalent noise power 
&w 0- £[ n ] i s tn e contribution of the residual 
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PN sequence in the received signal. From ( fl2l ). it can be seen 
that the linear convolution of the channel and data becomes a 
circular one after OLA process. Therefore, after fast Fourier 
transform (FFT), the received frequency domain TDS-OFDM 
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signal has finally a similar representation as CP-OFDM: 

k]=-=^yV[n]e-^ nk =H[i, k]X[i, k]-W'[i, k] (13) 

v n— 

where W is AWGN with the same variance as w'. The 
TDS-OFDM signal can thus be easily equalized by a one-tap 
equalizer: 

Y[i,k] 



Z[i,k] 



H[i,k]' 



(14) 



III. PN-based Channel Estimation 



As specified in |4|, the z^-length PN sequence in the GI is 
composed of an Apn -length PN sequence, more specifically 
a maximum-length sequence (m-sequence), as well as its pre- 
and post-circular extensions. Since any shift of an m-sequence 
is itself an m-sequence, the GI can also be treated as another 
ApN-length m-sequence with its CP. If this CP is longer than 
the length of CIR, the Apn -length PN sequence is ISI-free. 
Using the ISI-free PN sequence, a least square (LS) channel 
estimation is made for the i th OFDM symbol: 

Hi[i,k] = ^r = H[i,k] + ^^ ) 0<k<N PN , (15) 

where P and S are obtained through Ap^ -point FFT applied 
on the transmitted and received ISI-free PN sequences, respec- 
tively. The CIR estimate is 



h[i,l] = 



1 
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PX 



Npn-I 
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k=0 
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(16) 



Consequently, the A-length CFR estimation Hi is obtained 
by applying A-point FFT on h[i, I]. The mean square error of 
Hi is: 
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(17) 



As mentioned before, the channel estimation error results in 
interference to the OFDM data symbols. From the analysis 
given in the appendix, the power of the interference is com- 
puted as: 
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The part in the bracket of the above equation is only de- 
termined by the PN sequence and is constant for a given 
sequence. Thus, the power of the interference is determined 
by the length of channel L and the MSE of CIR estimation. 
Hence, the worse the channel estimation, the more interfer- 
ence will be introduced to the OFDM data symbols, which 
motivates us to propose additional processing to improve the 
channel estimation results. 

IV. Data- aided Channel Estimation 
A. Instantaneous Data-aided Channel Estimate 

In contrast to the classical turbo channel estimation al- 
gorithm like |fT3ll , the new proposed method excludes the 
deinterleaving, channel decoding and interleaving processings 
from the feedback loop as shown in Fig.Q] In other words, the 
soft data symbols used for data-aided channel estimation are 
rebuilt using the likelihood information from the demapper. 

The soft-output demapper demodulates the complex data 
symbols into Log-likelihood ratio (LLR) of bits 1 16 1. The LLR 
k] corresponding to the Z th bit of the [i, k)th equalized 
data symbol Z[i, k] is defined as: 



Xi[i,k] = log 



P(b[i,k,l] = l\Z[i,k}) 



P(b[i,k,l] =Q\Z[i,k]Y 



(19) 



where the P(b[i,k,l] = l\Z[i,k]) is the conditional proba- 
bility of the I th bit equal to 1 given Z[i,k]. The sign of 
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the LLR decides the corresponding bit equal to 1 or 0, 
and its absolute value gives the reliability of the decision. 
Based on the LLR, the probabilities of a bit equal to 1 and 

are P (&[*, k, 1} = 1) = and P (&[», k, 1} = 0) - 

1 — P(b[i,k,l] = 1), respectively. Then they are used as a 
priori probabilities to estimate the data symbols. First, the 
probability that the transmitted symbol X[i,k] is equal to a 
specific constellation point ctj is computed as the product of 
the probabilities of all the bits belonging to this constellation: 



log 2 fl 



P(X[i,k]= aj )= f[ P(b[i,k,l] = K l (a j )), 



(20) 



where ^ is the set of the constellation points of a given modu- 
lation scheme, [i is the modulation order and Ki(aj) € {0, 1} 
is the value of the I th bit of the constellation point otj. The 
soft data symbol is an expected value taking the a priori 
probabilities (|20| > into account: 

X[i,k]= a r P(X[i,k] = a j ). (21) 

Note that, being different from the classical decision feed- 
back channel estimation methods such as [10] and IH31 . no 
decision is made here in order to prevent error propagation. 
Directly using the soft data symbols, an instantaneous channel 
estimate is obtained over all active subcarriers by: 

H 2 [i,k] = — I— X[i,k]*Y[i,k] 



X[i,k]*X[i,k] 



>-H[i,k]- 



1 



V x 



■X[i,k]*W"[i,k](.22) 



where is the noise and interference component with 

a variance ct^„ = a^, + a], r]^[i,k} — \X[i,k]\ 2 is the 
power of the (i, fc)th soft data symbol which is used as power 
normalization factor. Note that, the power of the constellation 
is a constant and known value in the uniform power constel- 
lation cases such as BPSK or QPSK. Hence, the r)%[i, k] can 
be approximated by the power of the constellation r\ a in order 
to reduce computational complexity. If the data symbols are 
perfectly rebuilt, e.g. SNR is high, namely X[i,k] = X[i,k], 
turns to an LS estimator: 



H 2 [i,k]=H[i,k] 



W"[i,k] 
y/vx[i, k] 



(23) 



Unlike the turbo channel estimation, the proposed method 
does not include any error correction before rebuilding the 
data. Therefore the rebuilt soft data symbols, i.e. X's, are 
affected by the noise as well as the channel fades. It is quite 
possible that the instantaneous channel estimate ( 1221 is inac- 
curate and even erroneous for some subcarriers. Fortunately, 
as both the time delay spread and Doppler spectrum of the 
channel are limited, the CFR is highly correlated, i.e. almost 
identical within the coherence bandwidth and coherence time 
ifTTl . Moreover, the coherence bandwidth spreads over several 
adjacent subcarriers in the OFDM system with large FFT size, 
and the coherence time is normally longer than several OFDM 
symbol durations in low and medium velocity cases. This 
enables us to refine the data-aided channel estimation using the 



correlation of the CFR. Several channel estimate refinement 
approaches are proposed in the following sections. According 
to the time-frequency range that data-aided approaches work, 
they are catalogued into 1-D (frequency domain) and 2-D 
(time-frequency domain) ones. 



B. 1-D Refinements 

The proposed 1-D refinement approaches process all the 
instantaneous data-aided channel estimates within one OFDM 
symbol to get an improved estimation result at each time. Only 
the frequency domain correlation property of the channel is 
exploited. As the 1-D refinement approaches are carried out 
in a (OFDM) symbol-by-symbol fashion, they can track a fast 
variation of the channel and require minimum storage capacity. 
The OFDM symbol index i is omitted in this section for the 
sake of notation simplicity. 

1 ) Moving Average: Since the channel frequency response 
is almost identical within coherence bandwidth [17|, the most 
straightforward way to improve the data-aided estimation is to 
perform a moving average over the instantaneous channel esti- 
mates with a length less or equal to the coherence bandwidth. 
More specifically, the instantaneous channel estimates within 
a particular range, namely within a "window", are averaged 
to get a more precise estimate for the central position of the 
window. After the window sliding over all active subcarriers, a 
refined channel estimation is obtained. Concretely, the moving- 
averaged CFR for the fc th subcarrier is expressed as: 



^[ fc ] = Ji E 



me&k 



H[k]y^X[m]*X[m) 
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X[m]*W"[m] 



H 



mee k 



X[m]*W"[m 



Vxl m \ 
< k < N - 1, (24) 



where Gh 



{fc-L^ij.A-L^J + i, 



is the set of subcarrier indices within the moving average 
window with the k th subcarrier its central frequency, M is the 
length of the moving average window and W" is the noise 
and interference component. The moving average length M is 
chosen less than or equal to the coherence bandwidth. It can be 
either empirically pre-selected according to the "worst case" 
or be adaptively chosen by computing the coherence band- 
width using the estimated CIR from initial PN based channel 
estimation. The variance of the data-aided CFR estimate for a 
specific subcarrier is: 

<T 2 R2 {t,k]=E{\H{i,k}-H 2 {i,k}\ 2 } 

= E {|— ]T ^ l Mw"[iM \ 2 \ 



ri x [i,m\ 



\X[i,m}\ 2 _°lr» 
M- ^ r n ji [i,m} 2 M 2 



A/f2 



E 



w k 



(25) 
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The MSE of the data-aided CFR estimate for the i th OFDM 
symbol can be estimated by averaging (l25T > over N subcarriers: 

JV-1 
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Specifically, given the uniform power constellation with a 
power of r/ a , 
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H 2 ,XD 
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(M Vc 
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i27v"E E 
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Mr/a 



fjjt, (27) 



where 77^. = i ^fcLo 1 l-^M 1 2 * s ^e mean power of the rebuilt 
data symbols. 

2) Wiener Filtering: The proposed Wiener filtering based 
refinement is depicted in FigH] It consists of two steps: 
averaging and interpolation. 

A bunch of subcarriers are first selected as "virtual pilots". 
The virtual pilots are equally spaced and the spacing Lf (in 
terms of subcarriers) is selected substantially fulfilling the 
sampling theorem ifTHl . More concretely, with a 2x oversam- 
pling ratio, Lf is determined so that LfL/N < 1/4. Denote 
S the set of virtual pilots indices, i.e. S = {fc p |fc p = pLf, < 
p < L77-J — 1}- The cardinality of EE (the quantity of virtual 
pilots) is Kf. 

The instantaneous channel estimates within the coherence 
bandwidth around a virtual pilot k p £ 5 are averaged to get a 
more accurate channel estimate H 2 [kp] . Its expression can be 
obtained by setting k = k p in d24l) . Repeating the averaging 
process in all virtual pilot positions, we can get Kf accurate 
CFR estimation samples H2 [kp] , V/c p £ S. 

To obtain the TV-length CFR estimation H2, a Wiener 
filtering lTT8l based interpolation is performed on these CFR 
estimation samples: 



H 2 [k] 



ojf[k,kp]H 2 [kp], Q<k<N-l, (28) 



where cJ/[fe, k p ]'s are the Wiener filter coefficients com- 
puted to minimize the estimation MSE E{| H(k) — 
H2(k) | 2 }. Rearrange the samples in vector form H2 



12 = 



[H 2 [ko], H 2 [ki], . . . , H 2 [kKf-i]\ ■ Accordingly, the output 
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Fig. 5. 2-D averaging and interpolating. The averaging is performed in the 
2-D averaging region (shaded area). The interpolation is first carried out in 
the frequency domain ((T)) and consequently in the time domain (©). 



vector is H 2 = [H 2 [0] , H 2 [1] , . . . , H 2 [N- 1]] and the vector 

form of CFR is H 2 = [H[0], H[l], . . . , H[N - 1]] T . The 
coefficients of Wiener filter utf are expressed as: 



(29) 



where a?/ is the Kf x N coefficient matrix of Wiener filter 
with the (m, n)th element u>f[m,n], = -. {II ^ 1 1 , } = 
(R?' +a 2 jj\) is the Kf x Kf autocorrelation matrix of CFR 
estimation samples on pilot positions obtained in d24k is 
the variance of the estimation error H 2 derived in ( |26] |. I is the 
identity matrix. 8 f = E{H 2 H} = R^ 2) is the Kf x N cross- 
correlation matrix of CFR estimates H 2 and the real CFR H. 
R^ 1 ' and R^ 2) are the autocorrelation matrix of the real CFR 
in the frequency domain, and its components can be computed 
by (0. Note that the coefficients can be computed in either 
non-adaptive lfT8l or adaptive manners [19|. The MSE of the 
CFR estimate after frequency domain filtering is: 



£ H 2 ,1D 



N V f 



Rf )T *?Rf> 



(30) 



where Tr(-) is the trace operation, 
autocorrelation matrix of CFR. 



R 



(3) 



is the N x N 



C. 2-D Refinements 

As presented in above, 1-D refinement approaches only 
exploit the frequency domain correlation property of the 
channel. More processing gain can be acquired if the time 
domain correlation property is also taken into account. 

1 ) 2-D Moving Average: Expanding the average range in 
(T24l to 2-D, namely to several consecutive OFDM symbols, 
more instantaneous channel estimates can be involved in the 
moving average process. With a proper selection of the averag- 
ing length in the time domain with respect to the time variation 
of the channel, the 2-D moving averaged results are expected 
to be better than the 1-D counterpart. More specifically, the 
subcarrier indices within the 2-D averaging region with the 
fc th subcarrier of the i th OFDM symbol its central position 
is given by the set 0^ = {{p, q)\i - [^^-\ < p < 
i+ [^Ik- < q < k+ L^r 1 ]}- The 2-D 
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moving-averaged CFR estimates are written as: 

L ' J M t Mf ^ r) ± \pA 

The variance of the estimation error in OP) is then: 

4 2 [z,fc] = E{|#M]-iJ 2 M]| 2 } 

l\M t M fp j£ ikVjt [i,k] I 
a 2 „ 1 

The MSE of the 2-D moving averaged CFR is: 

JV-l 
fc=0 

= w" V V 1 (33) 

2) 2-D Wiener Filtering: The 2-D estimate refinement and 
interpolation is carried out over a number of OFDM symbols 
within an interpolation block as shown in Fig. [5] The block 
size B, i.e. the number of OFDM symbols in the block, is 
chosen according to the constraints of latency, storage and 
complexity. 

For low-complexity consideration, some time-frequency in- 
dices (i p ,k q ) £ f2 are pre-selected as virtual pilots. Lf and 
Lt are the virtual pilot spacings in the frequency and time 
domains, while Kf = and K t = are the number 

of virtual pilots within one OFDM symbol and the number 
of OFDM symbols including virtual pilots, respectively. The 
selection of virtual pilots in the frequency domain is discussed 
in the previous section. The selection in the time domain 
should also satisfy the sampling theorem. In a similar manner 
as in the frequency domain, the spacing between virtual 
symbols in the time domain L t (in terms of OFDM symbols) 
is selected so that L t Tbfd < 1/4. Hence, the set of indices of 
virtual pilots is ft = {(i P) k q )\i p = pL t , < p < K t ~ l;k q — 
qL f ,0<q<K f -l}. 

Then a refined CFR estimate for the virtual pilot located 
at (ip,k p ) is computed by averaging all the available CFR 
estimates within the coherence region ^ . The averaged 
estimate can be obtained from (|3TT > by setting i = i p and 
k = k p . Repeating this process over all virtual pilots in 
the interpolation block, K t ■ Kf refined CFR estimates are 
acquired. With these more reliable estimates, the overall CFR 
estimate can be obtained via interpolation. 

The 2-D Wiener filtering outperforms other interpolation 
techniques, e.g. linear interpolation, FFT-based interpolation 
etc. [20|. The use of two concatenated 1-D Wiener filters, 
i.e. one in the frequency domain and the other in the time 
domain, significantly reduces the computational complexity 
with negligible performance degradation compared to the 



optimum 2-D Wiener filter lfl"8l . The 2xl-D Wiener filtering 
is represented as: 

H 2 [i,k] = )^oj t [i, k, ip] Uf [k, i p , k p ] H 2 [i P , k p ] , (34) 

» ' 

time domain frequency domain 

where w t 's and w/'s are the coefficients of the 1-D Wiener 
filter in the time and frequency domains, respectively. The 
coefficients of the frequency domain Wiener filter are the 
same as those given in ( 129) . Let us define the filtered channel 
estimate matrix H2 of size Kf x K t whose (m,n)th element 

is H2[m, n). The output of the frequency domain 1-D Wiener 

„ f _ 

filter is denoted as an N x K t matrix H 2 = u>jH2, where u>/ 
is given in (1291 . The corresponding MSE of the estimation 
after frequency domain Wiener filtering is shown in ( T30b . 

The coefficients of the time domain Wiener filter are com- 
puted as: 

uf t = (35) 

where # t = EjH^H^*} = (r| 1} + ct 2 - / I) is the K t x K t 

autocorrelation matrix of the filtered CFR estimates H 2 , a^ f 
is the variance of the estimation error derived in d30l ). 6 t = 
E{H 2 fe H*} = R 4 (2) is the K t X B cross-correlation matrix of 

CFR estimates for the fc th subcarrier with different time H 2 k 
and the real CFR for a particular subcarrier with different time 
H t . Finally, R^ 1 ' and r[ 2) are the autocorrelation matrix of the 
real CFR in the time domain, and its components are computed 
by ©. 

The MSE of the CFR estimate after time domain filtering 

is: 

£42fl = i Tr (Rf'-^$-Rf"). (36) 

According to the Jakes' model, the coefficients of the time 
domain Wiener filter are only determined by the Doppler 
frequency. Note that, since the time domain interpolation 
length is not very large, the computation of the coefficients of 
time domain Wiener filter is not prohibitive. It should be also 
noted that the same coefficients are used for all subcarriers, 
which indicates that the storage spent for these coefficients is 
negligible compared with the size of interleaves 

Concerning the frequency domain Wiener filter, the coef- 
ficients can be pre-computed according to the uniform delay 
power spectrum assumption which represents the "worst case" 
of the mobile channel lfT8l . Given a virtual pilot pattern, 
the values of the coefficients only depend on the maximum 
delay spread of the channel which can be measured in reality. 
Therefore, the coefficients can be pre-computed with some 
typical channel situations. Even more complexity reduction 
can be achieved by limiting the Wiener filtering on several 
neighboring pilots [20). 

D. MMSE Combination 

When both PN-based and data-aided channel estimates are 
obtained, a linear combination is proposed to get a final CFR 
estimate: 

H = (3H 1 + (1-/3)H 2 , (37) 



TABLE I 

Complexity and Storage Comparison 



Method 


Complexity 


Storage 


PN based method 


• lno- lA 


jy rnmnlpv svmhnU 


Proposed 
method 


Moving average 


O(N) 


One OFDM symbol for 1-D methods, B 
OFDM symbols for 2-D methods. 


Wiener filtering 


O(N) to U(N ) 


MUCK03 [2| 


0{v ■ logv) when M -C v, 0(u 2 ) otherwise 


M OFDM symbols a 


YANG08 [TT] 


0(K ■ log/C) b 


One OFDM symbol and one GI 


WANG05 Q] 


0(N ■ log N) 


Two OFDM symbols 


TANG07 (9) 


0({N + u)-log(N + u)) 


One OFDM symbol and one GI 


STEENDAM07 (3) 


0(N 3 ) 


One OFDM symbol 


ZHAO08 Q3| 


(O(N) + complexity of decoding ) for iteration, 
0(N 2 ) or 0(N 3 ) for final stage 


Storage should be greater than the interleav- 
ing depth. c 



a The parameter M is the length of the averaging window. It is set to 20 to 40 OFDM symbols for BPSK and QPSK, 40 to 72 

OFDM symbols for 16QAM and 120 to 240 for 64QAM. 
b K, is the size of FFT/IFFT used to build the received PN sequence and make channel estimation. It can be selected from v + L to 

N or even larger. K. is set to 2048 in II 11 . 
c As far as the DTMB system is concerned, the interleaving depth is 170 or 510 OFDM symbols. 



where _ff 2 is a generic expression of the data-aided channel 
estimations which can be obtained from (1241 . ( f28l l. (l3~TT l or 
(l34l using different techniques. The combination is carried out 
only for the active subcarriers while the CFR estimates for null 
subcarriers are kept as the results obtained from the PN-based 
one. The optimum weight values (3 opt can be obtained using 
the MMSE criteria: 

P opt =argmm{E{|tf- J ff| 2 }} 

= argmmj/^ + Cl-/}) 2 ^}. (38) 

The above equation uses the fact that H\ and H2 are obtained 
from different sources and thus uncorrelated. Since the MSE 
function is convex, the f3 opt is obtained by setting the deriva- 
tive with respect to (3 equal to zero and the solution is: 

Po P t = — ^ , (39) 

£ Ht + £ H 2 

where Ef, is MSE of the PN-based estimation obtained 

ill 

from ( U7j >. while is the MSE of the data-aided channel 
estimation and is computed from (l26l l. (l30l >. d32l or (|36| > 
depending on different methods. 

The proposed method can be carried out in an iterative 
manner. More precisely, if it has not achieved the preset 
maximum iteration times yet, the combined CFR estimates 
(l37l i are used as the initial channel estimates of the next 
channel estimation iteration for PN subtraction and symbol 
equalization. The equalized data symbols using the updated 
channel estimates are more accurate than previous iterations. 
Hence, more reliable data-aided channel estimates can be 
acquired in the new iteration. If the the maximum iteration 
times are achieved, the combined channel estimates output as 
the final channel estimates for the data equalization and FEC 
decoding. 

V. Complexity Analysis 

In this section, we analyze the additional computational 
complexity introduced by the proposed data-aided channel 



estimation method. The complexity is evaluated in terms of 
required real multiplications and real additions for each OFDM 
symbol per iteration. In this paper one complex multiplication 
is counted as 2 real additions and 4 real multiplications al- 
though there exist some smarter ways requiring fewer multipli- 
cation times. The data rebuilding process including ( fl9l ), ( f20b 
and (1211) needs (N ■ /i log 2 fi+N-2 log 2 fi) multiplications and 
N- 2 log 2 /i additions. The instantaneous data-aided estimation 
d22l needs 87V multiplications and 3iV additions. 

The computational complexities required by the estimate 
refinement vary with respect to different strategies. For the 1- 
D moving average method, the averaging for each subcarrier 
with an averaging length of M requires 2{M — 1) additions 
and 2 multiplications neglecting the edge effect. Hence, the 
moving averaging over N subcarriers requires 2 (M — 1)N 
additions and 2N multiplications. The combination of the two 
channel estimates requires 2N additions and AN multiplica- 
tions. The overall basic operations are + 2) log 2 /i + 14) iV 
multiplications and (2 log 2 /1 + 2M + 3)N additions for each 
iteration. Given M <C N, the additional complexity of the 1-D 
moving average based channel estimation is O(N). 

The 2-D moving average method, on the other hand, has 
generally the same computational complexity except that the 
averaging process for each subcarrier requires 2(M t Mf — 1)N 
additions and 2N multiplications and the overall required 
additions turns to (21og 2 /i + 2M t M f + 3)N. Given the 
condition M t M f < N fulfilled, the complexity is still <D(N). 

As far as the Wiener filtering based approaches are con- 
cerned, for the K virtual pilot in each OFDM symbol, the 
averaging of the instantaneous estimates in the coherence 
bandwidth requires 2{M — \)K real additions and 2K real 
multiplications. As the coefficients of the Wiener filter is 
real [18], the Wiener filtering needs 2KN multiplications and 
2KN additions. Taking the symbol rebuilding and combi- 
nation into account, the overall operations for 1-D Wiener 
filtering based method is (2 log 2 fi + 5 + 2K)N + 2(M - l)K 
additions and ((fi + 2) log 2 fi + 12 + 2K)N + 2K multiplica- 
tions, which is between O(N) and 0(N 2 ). 

For the 2-D Wiener filtering case, in each block consisting 
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SNR (dB) 

Fig. 6. MSE of the proposed data-aided channel estimation using 1 -D Moving 
Average and 1-D Wiener filtering with QPSK, in the TU-6 channel and SFN 
channels with velocity of 30 km/h. 

of B OFDM symbols, the averaging and 1-D frequency 
domain Wiener filtering is repeated K t times which counts 
2(N+M t M f -l)KfK t additions and 2(N+l)K f K t multipli- 
cations. Then, the time domain Wiener filtering needs 2K t BN 
multiplications and additions. Therefore, the overall operations 
per OFDM symbol per iteration for 2-D Wiener filtering based 
method is (((J.+2) log 2 fi + l2 + 2^ + 2K t ) N multipli- 
cations and (2 log 2 n + 5 + 2^ + 2K t )N + 2 iMlMl^SlEl 
additions, which is, similar to the 1-D Wiener filtering case, 
between O(N) and 0(N 2 ). Even lower complexity can be 
achieved by reducing the size of the Wiener filter. 

The computational complexities of proposed method and 
several typical existing methods 0, (3), Q, (9), flT], £13] 
are shown in Table U From the comparison, we find that 
the method MUCK03 achieves least computational complexity 
but requires a large amount of storage. The proposed method 
using moving average needs less complexity than rest of 
other methods. In the meantime, the complexity and storage 
costs of proposed method using Wiener filtering are close to 
the methods WANG05, TANG07 and YANG08. The methods 
STEENDAM07 and ZHAO08 spend more complexity than other 
techniques. 

VI. Simulation Results 

A. Simulation Settings 

Simulation parameters are chosen from the DTMB sys- 
tem [4|. Baseband signal bandwidth is 7.56 MHz. FFT size 
is 3780 which results in a subcarrier spacing of 2 kHz. All 
subcarriers are active. The length of the GI is set to v — 420. 
The power of the PN sequence is twice as much as data sym- 
bols. The COST207 Typical Urban (TU-6) channel model ED 
is employed in the simulation. The proposed algorithm is 
also evaluated in the single frequency network (SFN) scenario 
which is a spectrum efficient solution widely used in broadcast 
networks. In this case, the same signal is sent from two 
different transmitters at the same time on the same carrier 
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Fig. 7. MSE of the 1-D Moving average method before combination with 
different averaging lengths. 

frequency. The signal from these two transmitters experiences 
independent fading. The equivalent CIR of the SFN channel 
is the combination of the CIR's of two independent TU-6 
channels. The propagation distance difference between the two 
signals causes time delay and power attenuation on the second 
CIR [14|. In this paper, the distance difference is set to 7 
km corresponding to 23.33 [is time delay which makes the 
overall length of the CIR longer than that of the CP of the 
PN. The power attenuation is set to 10 dB. Following the 
computation in [17], when the frequency correlation function 
is 0.9, the coherence bandwidths of the TU-6 channel and 
the SFN channel are 18.8 kHz and 2.94 kHz, respectively. 
Given the 2 kHz subcarrier spacing in the DTMB system, the 
frequency domain averaging length M is accordingly set to 9 
and 3 for the TU-6 and SFN channels, respectively. The virtual 
pilot spacings in the frequency and time domains are set equal 
to the averaging lengths in each domain, respectively. 

B. MSE Performance 

FigH] presents the MSE performance of the proposed al- 
gorithm with QPSK, using 1-D moving average and 1-D 
Wiener filtering in the TU-6 and SFN channels, respectively. 
It can be observed that the accuracy of the channel esti- 
mation is progressively improved. The iterative estimation 
process converges very fast. A significant improvement can 
be obtained after only two iterations. More specifically, in 
the TU-6 channel, when the 1-D moving average is used, 
the proposed data-aided channel estimation method acquires 
about 4.1 dB gain in terms of required SNR to achieve 
an MSE level of 1 x 10~ 2 compared with the PN-based 
one. If the more powerful 1-D Wiener filtering technique is 
adopted, the gain increases to 5.1 dB. The higher the SNR, 
the more efficient the Wiener filtering. Moreover, it should be 
noted that the performance of using Wiener filtering after one 
iteration is as good as the best performance that using moving 
average technique can achieve. This provides a less processing 
delay (less iterations) trade-off with higher computational 
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SNR (dB) 

Fig. 8. MSE performance of the proposed method after two iterations using 
2-D moving average and Wiener filtering methods. Mt = 2, Mr = 9 for the 
TU-6 channel and Mt = 2, Mf = 3 for the SFN channel. Velocity is set to 
6km/h for both TU-6 and SFN channels. 



complexity for each iteration. On the other hand, in the SFN 
channel, since the delay spread is much longer than the CP of 
the PN sequence, there is a strong ISI on the PN sequence. 
This can be observed from the fact that the performance of 
the PN-based channel estimation is seriously degraded and 
appears an estimation error floor at high SNR. Furthermore, 
as the averaging length is much shorter in the SFN case, the 
noise mitigation ability is limited in the data-aided channel 
estimation. Even though in such a harsh channel, the data- 
aided method however offers 6.9 dB and 8.1 dB gain in terms 
of required SNR to achieve MSE of 5 x 10~ 2 compared to 
the PN based method, when using the 1-D moving average 
and Wiener filtering, respectively. Moreover, the estimation 
performance is improved approximately ten times at high SNR 
region. 

FigGl depicts the impact of the averaging length to the 
performance of 1-D moving average based method. At a low 
SNR, e.g. 10 dB, the MSE is monotonically decreasing with 
the increase of the averaging length. However, at a higher 
SNR, e.g. at 30 dB, the MSE of the averaged estimation 
results is degraded with a long averaging length. Hence, a 
moderate averaging length, say M — 7 ~ 9, is a proper choice. 
Considering that the data-aided channel estimation is more 
crucial in lower SNR region, it is better to bias our selection 
to a greater length. Therefore, M = 9 is a good trade-off that 
suits all noise level for the TU-6 channel, which coincides with 
the selection according to the coherence bandwidth. Similar 
conclusion can be drawn in the SFN channel. M = 3 is the 
best trade-off of the averaging length in the SFN channel. 

FigJS] shows the MSE of proposed 2-D method after two 
iterations. As the 1-D method is proved to be effective with 
QPSK, the 2-D method is only evaluated with higher order 
constellations, namely 16QAM and 64QAM. The averaging 
length in the time domain is selected to M t = 2 which 
is a very conservative setting that only requires the channel 
keeping similar within two consecutive OFDM symbols. Even 
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Fig. 9. MSE performance of the proposed method using 2-D Wiener filtering 
method after two iterations with different velocities. Mt = 2, Mf = 9 for 
the TU-6 channel and M t = 2, M f = 3 for the SFN channel. 

with such a short averaging length in the time domain, the 2- 
D method can still provide effective estimation. For instance, 
the proposed method using 2-D Wiener filtering acquires 8.4 
dB gain over the PN-based method in terms of the required 
SNR to achieve an MSE level of 1 x 10~ 3 with 16QAM in the 
TU-6 channel. In the SFN channel, the improvement is also 
significant. The estimation error floor is reduced from 2 x 10 -2 
to 3 x 10~ 3 with 16QAM at a SNR of 30 dB. Given longer 
time domain averaging length, more gain can be expected in 
a static channel. 

Fig |9] presents the changing of the MSE performance with 
different velocities of the receiver in the TU-6 and SFN 
channels. With the increase of the velocity, the variation of the 
channel among consecutive OFDM symbols is more and more 
notable, which limits the accuracy in the averaging results 
in (ED- In low SNR region (e.g. SNR less than 10 dB), the 
channel variation is a less significant influence that affects the 
estimation performance compared to the noise. Whereas in a 
higher SNR region where the noise is no longer the dominant 
factor, the channel variation compromises the performance of 
the proposed method. However, it should be noted that the 
proposed method can still provide adequate improvement. For 
instance, in the TU-6 channel with a velocity of 30 km/h, the 
proposed method obtains 8.3 dB gain with 16QAM in terms 
of required SNR to achieve MSE of 1 x 10~ 3 . While in SFN 
channel, as the ISI is the dominant factor that affects the MSE 
results. Therefore, the impact of the channel variation is not 
significant. 

C. BER Performance 

In this subsection, we present the bit error rate (BER) per- 
formance of the DTMB system using the proposed data-aided 
channel estimation method. Both the convolutional interleaver, 
BCH code and LDPC code are included in the simulation 
in order to give a system level evaluation. The convolutional 
interleaver is set to (52, 240) which corresponds to a time 
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Fig. 10. BER of the DTMB system with QPSK and LDPC (R=0.8) using 
different channel estimation methods in the TU-6 channel with velocity of 30 
km/h. Iteration time is set to 2 for the proposed data-aided channel estimation 
in both 1-D Moving Average (MA) and 1-D Wiener Filtering (WF) cases. 



Fig. 12. BER of DTMB system with 64QAM and LDPC (R=0.6) using 
different channel estimation methods in the TU-6 channel with velocity of 
30km/h. M t = 2, Mf = 9 for Moving Average (MA) and Wiener Filtering 
(WF) cases. 
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Fig. 11. BER of the DTMB system with QPSK and LDPC (R=0.8) using 
different channel estimation methods in the SFN channel with velocity of 30 
km/h. Iteration time is set to 2 for the proposed data-aided channel estimation 
in both 1-D Moving Average (MA) and 1-D Wiener Filtering (WF) cases. 
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Fig. 13. BER of DTMB system with 64QAM and LDPC (R=0.6) using 
different channel estimation methods in the SFN channel with velocity of 30 
km/h. Mt = 2, Mt = 3 for both moving average and Wiener filtering cases. 



delay of 170 OFDM symbols. The mobile speed is set to 
30 km/h which introduces a maximum Doppler frequency of 
13.89 Hz, given the carrier frequency 500 MHz. The coherence 
time is accordingly 0.0305 s and equivalent to the duration 
of 61 OFDM symbols (without GI). The proposed method is 
compared with the PN-based method as well as three important 
methods in the literatures, namely MUCK03, TANG07 and 
YANG08. Since the MUCK03 method is the least complex one 
in the literature and the TANG07 and YANG08 methods have 
comparative complexity as the proposed one while providing 
superior performance than its predecessors such as WANG05, 
it is interesting to take them as reference methods. 

Figd0]and Fig[TT]present the BER performance with QPSK 
in the TU-6 and SFN channels, respectively, while Figfl2l 
and Fig [JJ] depict the BER with the 64QAM in these two 



channels. Note that the method TANG07 does not adapt to 
the SFN channel. Therefore, its performance is not shown in 
that case. As the method MUCK03 is affected by the channel 
variations, its performance is presented with low mobility (6 
km/h) and medium mobility (30 km/h) cases. To quantify the 
performance of each method, we observe the required SNR 
to achieve a BER level of 5 x 10~ 5 . The gains obtained by 
the proposed method are listed in Table HU There are some 
interesting observations to be highlighted: 

1 . The proposed method outperforms the reference methods. 
The proposed method achieves better BER performance than 
the reference methods in almost all cases considered in the 
simulations. More specifically, it acquires more than 0.5 
dB, 0.3 dB and 1.1 dB gains over the PN-based, MUCK03 
and TANG07 methods, respectively, with QPSK in the TU- 
6 channel. When used with 64QAM, the improvements are 
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TABLE II 

Summary of the gains obtained using the data- aided channel estimation methods. 



Channel 


Modulation 


Method 


Gain over 
PN-base 
memou \cld ) 


Gain over 
method 

1 AINOU / \CLD ) 


Gain over 
method 

IVIUO(\UO \CLD ) 


Gain over 
method 

YrtlNuUO \dD ) 


Gap away from 
perfect estima- 
lion (Qd j 


TU-6 


QPSK 


Moving average 


0.5 


1.1 


0.5 [0.3] 


-0.05 


0.5 


Wiener filtering 


0.6 


1.2 


0.6 [0.4] 





0.4 


64QAM 


Moving average 


0.3 


0.5 


- [> 3.0] 


-0.05 


0.3 


Wiener filtering 


0.4 


0.6 


L - [> 3.1] 





0.2 


SFN 


QPSK 


Moving average 


1.7 




> 4.3 [-0.1] 


1.38 


1.3 


Wiener filtering 


1.8 




> 4.4 [0.04] 


1.46 


1 2 


64QAM 


Moving average 


> 7.0 






> 2.0 


5.8 


Wiener filtering 


> 8.0 






> 3.0 


4.2 



The first value is obtained with M = 20 given a velocity of 30 km/h, while the value in the square brackets is obtained with 
M = 50 given a velocity of 6 km/h. 



accordingly more than 0.3 dB, 3 dB and 0.5 dB, respectively. 
In the more harsh SFN channel case, the improvements are 
more obvious. The proposed method provides 1.7 dB and 4.3 
dB gains over the PN-based method and MUCK03 method 
with QPSK, respectively. Compared with YANG08 method, 
the proposed method achieves similar performance in the 
TU-6 channel. Yet, we also notice that after the partial- 
decision-aided interference removal there still exists residual 
interference in the PN sequence. The impact of the residual 
interference is not serious when the channel delay spread 
is short (for instance in the TU-6 channel). However, it 
becomes more harmful in the SFN channel. The performance 
of YANG08 method significantly degrades in SFN cases. The 
proposed method outperforms YANG08 method more than 1.3 
dB and 2.0 dB with QPSK and 64QAM, respectively. In 
addition, the residual gap between the proposed method and 
the perfect channel estimation case turns out to be very small. 
Especially in the TU-6 channel, the proposed method is only 
0.2 dB to 0.5 dB away from the case using perfect channel 
estimation. 

2. The proposed method is robust in different channel condi- 
tions. It can steadily provide satisfactory performance in time- 
varying channels as well as the channels with long time delay 
spread. The method MUCK03 eliminates the ISI by averaging 
the known sequence over a large number of OFDM symbols 
and thus requires the channel to remain unchanged within the 
averaging window. It is consequently not robust in face of the 
time-varying channel condition. In addition, it is not efficient 
with higher order constellations such as 64QAM. In contrast, 
the proposed method performs channel estimation over only 
either one (1-D case) or two OFDM symbols (2-D case). 
Therefore, it has stronger immunity to the channel variations. 
Moreover, the SFN channel has extremely long time delay 
spread with strong echoes and is thus recognized as the most 
difficult situation for the channel estimation. For instance, the 
method TANG07 suffers from the interference from previous 
OFDM symbols and does not adapt to the SFN channel. The 
performance of YANG08 method significantly degrades in the 
SFN channel due to the interference as well. In contrast, the 
proposed method exploits the OFDM data symbols to refine 
the channel estimates obtained from the PN-based estimation, 
which can provide reliable channel estimation in presence of 



long ISI. 

3. Both the moving average and the Wiener filtering tech- 
niques provide satisfactory performance with different channel 
conditions. The performance difference between two meth- 
ods is small in the TU-6 channel and in the SFN channel 
with QPSK. That is to say it is convenient to use the low- 
complexity moving average based estimation in the common 
case. However, the Wiener filtering technique performs better 
with 64QAM in the harsh SFN channel where it requires more 
reliable channel estimation. Moreover, it has been shown in 
the MSE simulation that the Wiener filtering technique needs 
fewer iterations to achieve a satisfactory estimation results. In 
other words, it is suitable for the application that has higher 
restriction of processing time delay. 

Note that all the results are obtained using the powerful 
LDPC code. More gains can be expected when weaker channel 
codings are used. 

VII. Conclusion 

This paper proposes a novel efficient solution to the chal- 
lenging channel estimation task in the TDS-OFDM system. 
The presented algorithm takes basis on the classical PN-based 
estimation and processes estimation refinement by exploiting 
a low-complex data-aided approach. In contrast to a classical 
turbo channel estimation case, this new algorithm does not 
include the decoding and interleaving functions when rebuild- 
ing the data symbols in order to reduce the computational 
complexity. This enables the new proposed method to be 
integrated in systems which contain extremely long interleaver 
and sophisticated channel decoder like the DTMB system. The 
rebuilt soft data symbols are used as virtual pilot for data-aided 
channel estimation. In order to suppress the noise existing in 
the data-aided channel estimates, several techniques including 
moving average and Wiener filtering in 1-D and 2-D have been 
conducted. The cooperation of these refining techniques and 
the soft symbol feedback enables reliable data-aided channel 
estimation without intensive computation. Simulation results 
show that the new proposed data-aided channel estimation 
algorithm provides satisfactory performance even often very 
close to the perfect channel estimation case. Compared to 
already existing techniques, the proposed method performs 
very well in various transmission cases, including the 64QAM 
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constellation, long delay spread channels and mobility scenar- 
ios. Note that the use of the novel low-complex data-aided 
channel estimation algorithm presented in this paper is not 
limited in TDS-OFDM. In fact, it can be easily adopted by the 
traditional CP-OFDM systems where the pilot assisted channel 
estimate will serve as the initial estimate. 

Appendix 

Interference Resulted from Imperfect PN 
Removal 

Take ( fTOb and ( fTTT i. and suppose that the PN sequences are 
identical for all OFDM frames, the residual PN sequence in 
the received signal is written as: 

L-l 

f [n] = ^2 Ah i4 n ~ 0<n<v. (40) 

The interference for the k th subcarrier is thus: 



N-l 



i[k] = -L= E £['■ 

v n— 
^ L-l 

= -t^E^E^-^ 



(41) 



;=o 



n=0 



Denote c[n — Z]„ by c;[n], representing a circular shift of I bits 
of sequence c[n] towards right side. 

Using the US assumption, i.e. E[h*h r> ] = of <S[Z — /'], the 
power of the interference is: 



crf[fc] = E{/[fc]*7[fc]} 



L-l 



= ^E^ E c 'Kr c iH^ f 



(«2— rai)fe 



i=0 



™1 ,^2—0 



where o^ ft = E[|Aft, z | 2 ] is variance of the estimation error 
of the I th channel tap. Substituting ri\, n-i by q = — fi\ and 
n = m, it yields: 
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E^< 

9=1 



^ E ^Ah, 



i/-i 
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,n=0 
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(42) 
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